Diffraction effects on light-atomic ensemble quantum interface 
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We present a simple method to include the effects of diffraction into the description of a light- 
atomic ensemble quantum interface in the context of collective variables. Carrying out a scattering 
calculation we single out the purely geometrical effect and apply our method to the experimental 
relevant case of Gaussian shaped atomic samples stored in single beam optical dipole traps probed 
by a Gaussian beam. We derive simple scaling relations for the effect of the interaction geometry 
and compare our findings to the results from one dimensional models of hght propagation. 



I. INTRODUCTION 



Coupling collective variables of atomic ensembles to 
propagating modes of the electromagnetic field has been 
identified as an efficient tool to engineer the states of 
atoms and light at the quantum level. Several proposals 
for spin-squeezing, mapping of quantum states between 
light and atoms, i.e. quantum memory operations, cre- 
ation of macroscopic entanglement, and teleportation of 
atomic states have been published [ISIEIllEIEIllll- 
A number of these proposals have actually been verified 
experimentally using atomic samples stored in vapor cells 
as well as in magneto optical traps 0, 0i 

hj. The 

efficiency of the coupling is often discussed resorting to 
effective 1-D models for the light propagating through 
a homogeneous sample and the optical depth is found 
to be the essential parameter determining the coupling 
strength. This naturally suggests the use of cold and 
trapped atomic samples, where long coherence times for 
collective variables are possible, with a shape mimicking 
a 1-D string of atoms to maximize column density for 
a given number of atoms. However, while the effective 
1-D models work well for wide and nearly homogeneous 
samples the question of light diffraction from small, in- 
herently inhomogeneous samples and its impact on cou- 
pling efficiency immediately arises when cold and trapped 
atomic samples are used (l^. In a recent work [lj| as- 
pects of spatial inhomogeneity have been addressed via 
the introduction of new asymmetric collective variables. 

The purpose of the present paper is twofold. We 
present a simple and effective semi-analytic method to 
identify the spatial light mode the atomic sample cou- 
ples to and use this method to find optimum shapes for 
atomic samples. In addition we will derive simple scal- 
ing parameters describing the effect of the shape of the 
atomic sample which allow us to make a link to the estab- 
lished 1-D models and quantify the effect of diffraction 



on the coupling strength. Albeit we are ultimately in- 
terested in measuring and manipulating quantum states, 
we find it instructive to solve the underlying classical 
scattering problem first to single out the purely geomet- 
ric effects, and we try to avoid hiding useful practical 
information behind the much more intricate mathemati- 
cal formalism needed to solve the full quantum problem. 
Our approach delivers valuable information for designing 
experimental configurations, provides intuitive pictures, 
and may serve as a starting point for a more elaborate 
quantum theoretical treatment. 

The remainder of the paper is organized as follows. In 
Section 2 we review briefly the quantum description of 
two modes of light coupled coherently to atoms in terms 
of collective variables and introduce typical experimen- 
tal configurations. Using this framework we motivate the 
use of the coherently scattered power by the atomic sam- 
ple as the relevant parameter to optimize for the quan- 
tum coupling. Section 3 presents our model for calculat- 
ing the scattering mode and the scattering efficiency of 
the atomic sample for different geometries together with 
some remarks on the assumptions made and the range of 
validity of the model. In Section 4 we apply this model to 
atoms stored in a single beam optical dipole trap. In Sec- 
tion 5 we provide a link between our classical calculation 
and the results from effective 1-D models. Section 6 sum- 
marizes our results and points out possible extensions of 
our model. 



II. COLLECTIVE LIGHT-ATOM COUPLING 

The hermitian part of the interaction of a two-mode 
pulse of off-resonant light with an ensemble of atoms 
with two ground states residing in a container with cross- 
sectional area A and length L can be described by a 
Hamiltonian of the form (see e.g.jTi llslll^lr^ ): 
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The factor in front of the integral describing the strength 
of the coupling of the collective variables contains the 
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single-atom reactive response to radiation characterized 
by the cross-section a^, the detuning A in units of the 
natural linewidth T, and the cross-sectional area of the 
light mode governing the electric field strength per pho- 
ton occupying the mode. In eq.l^ collective operators 
for light S and atoms J are introduced. For light we 
can write S in terms of continuous mode creation and 
annihilation operators with boson commutation rules: 
Sx = i(ajai -I- 0+02), 5y = ^{afai - afa2),Sz = 
^(0^0.2 — a^fli). For simplicity of notation the time and 
spatial dependence of the operators has been suppressed. 
The mode index can specify for instance two fields prop- 
agating in the two legs of a Mach-Zehnder interferometer 
(Fig^) or left and right circular polarized modes of a 
single laser beam (Fig^), in which case the operator 
expressions describe the instantaneous Stokes vector of 
the light 13 ■ By construction the collective operators 
follow the usual commutator algebra for vector opera- 
tors. For atoms, restricting us for simplicity to the case 
of two ground-states only, we define in analogy: — 
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|1),(2|), J, = Sf-W i(|2),(2| - where the sum 

extends over the number of atoms Nat{z) residing in a 
thin slice of the interaction volume at position z .19j. A 
continuum limit is taken by shrinking the thickness 5z of 
the slices to zero. The three operators Ji are the macro- 
scopic Bloch-vector components for the 2-level system 
formed by the two ground states. Inherent in this descrip- 
tion are the assumptions that the transverse mode func- 
tion is a frozen variable, hence the 1-D integration over 
the Hamilton density in eq.(^, and that the atoms scat- 
ter light coherently only into the input modes. This is a 
good approximation for atomic ensembles and light fields 
transversally much wider than an optical wavelength and 
not too long samples. In addition, since the operator na- 
ture of atomic position is suppressed, the spatial density 
distribution of the atoms is not changed by the interac- 
tion with the radiation in this model. 

The dynamics generated by eq.lQ) at the level of ex- 
pectation values describes differential phase shifts of the 
light modes due to coherent forward scattering by the 
atoms in the two ground states and differential shifts of 
the atomic energy levels due to the intensity difference 
of the two light modes. This dynamics can be conve- 
niently visualized with the help of coupled Bloch spheres 
for light and atom variables. Of particular interest for 
practical quantum state engineering is the case of initial 
conditions with J™ and large and classical, i.e. un- 
correlated coherent states with expectation values on the 
order of half the total number of atoms and photons re- 
spectively. In this case the dynamical evolution of the 
system is driven entirely by the fluctuations of light and 
atomic variables. In the experimental setup sketched in 
Fig^ this corresponds to linearly polarized input light 
and all atoms in a balanced coherent superposition of the 
two ground states, such that there is a macroscopic mag- 
netic polarization in the plane orthogonal to the direction 



of propagation. Linearizing the spectrum of the opera- 
tors around their expectation values and integrating the 
dynamics, to leading order, input/output relations can 
be written for the fluctuations of the collective operators 
orthogonal to the coherent excitation [23, 0| : 
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From these relations we can calculate the variance of a 
set of measurements on identically prepared systems. A 
set of measurement records of 5°"* done over some time 
interval will show a variance: 
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Here the first term reflects the noise of the coherent in- 
put state of light, with riph the total number of probe 
photons, while the second term describes the excess fluc- 
tuations imprinted onto the light by the fluctuations of 
the atomic variable and amplified by the large coherent 
amplitude of the conjugate light variable. From the mea- 
surement record we can infer a value for the atomic vari- 
able with a confidence interval limited by the first term, 
i.e. with a signal to noise ratio of: 
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Projection of the atomic state by the destructive mea- 
surement on light, will leave the collective atomic variable 
in a state with fiuctuations reduced below the standard 
quantum limit if the signal to noise ratio is finite. In 
practical terms this means that repeated measurements 
on the same atomic system will show a covariance below 
the projection noise [23|. To the extent the measurement 
is nondestructive, i.e. after proper inclusion of sponta- 
neous emission [23|, this fact can be exploited to write, 
store, and retrieve quantum information to and from the 
atomic variable. Instead of measuring destructively the 
light pulse after a single pass through the atomic ensem- 
ble one can allow for multiple interaction together with 
appropriate switching of the coupling Hamiltonian be- 
tween consecutive passes. This way unconditional quan- 
tum state preparation is possible and the measurement 
serves merely as a verification of successful state prepa- 
ration I23 13 . In this more general setting the above 
defined signal to noise ratio plays the role of a coupling 
strength k^: 



(3) 




Figure 1: (color online) Mach-Zehnder and polarization interferometry setups for quantum coupling of collective variables of 
cold atomic samples to light. 



which can be conveniently expressed as the prod- 
uct of integrated spontaneous emission rate 
V = nph{<^o/A) (r/2A)^ and the optical depth or 
column density ao = Nat{<JQ/A), explaining the use 
of the optical depth as a figure of merit for collective 
variable light-atom coupling 24]. The interaction 
geometry, of course, stays the same, so diffraction of 
light discussed in the following has the same impact on 
the coupling strength in these schemes. 



pulse duration, in the presence of atoms is given by: 



After this brief review of the 1-D quantum model, we 
look again at the experimental configuration in Fig^ 
[23, seeking this time a description in purely classical 
terms and pointing out where we go beyond the 1-D 
model. A cloud of cold atoms residing in one arm of 
the interferometer is loaded into a far-off-resonant dipole 
trap created by a focused gaussian beam. Probe light 
with a wave- vector k = 27r/A enters the interferometer 
from the left and passes through the atomic sample. The 
light experiences a phase shift caused by the refractive 
index of the sample, which is determined by the atomic 
density and the distribution of population among the two 
different sub-levels. Given an atomic density distribution 
the light carries thus information about the level popu- 
lations after passage through the sample. In turn, the 
atomic energy levels are shifted differentially during the 
interaction and thus information about the light is de- 
posited in the coherence between the atomic levels. To 
assess the sensitivity of a measurement of the atomic pop- 
ulation difference we look at the signal by the balanced 
homodyne detector (with quantum efficiency e) on the 
right upon detection of a pulse of light (frequency lu) of 
duration r. The detector signal S'd, i.e. the photocurrent 
is in units of the elementary charge e integrated over the 
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Here Eprobe{Eref ) denotes the initial field in the probe 
(reference) arm. The integrals over the detector area 
describe the overlap of the transverse mode functions 
Msc, Mprobe, Mref in the detector plane and contain also 
the oscillatory dependence on the path length difference 
of the interferometer. We separate explicitly the field 
scattered by the atoms E^c from the probe field Eprote, 
since in general it will not be in the same spatial mode 
and have a different phase. We can adjust the path length 
difference such to make the second integral in eq. Q van- 
ish - this corresponds to the preparation of a large (S^) 
and a measurement of (Sy) with the detector in the lan- 
guage of the quantum description. We note that for dis- 
persive scattering in the far field, probe wave and scat- 
tered wave are nearly tt/2 out of phase, so the first term 
in eq.Q will take its maximum value whenever the sec- 
ond term vanishes. The detector signal stems then from 
the interference of the reference wave with the scattered 
wave only. The atomic contribution and consequently 
the sensitivity of the measurement will thus be a mono- 
tonic function of both the coherently scattered power by 
the atoms and of the geometric overlap of the wavefronts 
of scattered wave and reference wave. Carrying out the 
overlap integral assuming perfect matching of the refer- 
ence wave to the scattered wave, we write for the detector 
signal power S^: 
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The noise will be limited from below by the shot noise of 
the detected photons 12^. We write the corresponding 
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noise power iV^ as: 

= —2Pref 

Since we want to assess the sensitivity with reference to 
the fluctuations of the population difference in an uncor- 
related atomic sample we have to normalize our signal 
power to this variance, i.e. to a half of the total number 
of atoms. With this we arrive at an equivalent expression 
for the signal to noise ratio (SNR) or coupling strength 
as: 
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Here the coupling strength is expressed in terms of 
the - admittedly artificial quantity (see sectional - total 
scattered power per atom. We conclude that optimizing 
the sensitivity is equivalent to maximizing the scatter- 
ing efficiency and adapting the reference wavefront to 
the scattered wave. Trivially the scattered power can be 
increased by increasing the input power to the interfer- 
ometer. But here upper limits are set by saturation of the 
detectors and the requirement for low spontaneous emis- 
sion rate. We choose therefore to perform the search for 
optimal shapes of the atomic samples and of the probe 
beam in the next section under conditions of constant 
probe power. Given a maximum pulse energy the detec- 
tors can withstand, the integrated spontaneous emission 
rate can be set to a desired value by choosing the single 
atom scattering cross section via the detuning. In the 
next section we calculate the scattering efficiency and 
the scattered wavefront for selected geometries. 



III. SCATTERING MODEL 

In order to calculate the stationary diffracted field by 
the atomic density distribution wc model the sample as 
an ensemble of fixed (i.e. infinitely heavy) point scat- 
terers and make use of a Born approximation [2^ l28l | . 
Assuming fixed positions we neglect Doppler and recoil 
shifts. Since we expect coherent diffraction to occur 
mainly close to the forward direction and we are inter- 
ested in samples at ultralow temperatures, Doppler shifts 
will play a negligible role. On time scales short compared 
to the recoil time - the time needed for an atom to travel 
a distance of an optical wavelength at one recoil veloc- 
ity - we can assume that spatial correlations are frozen, 
i.e. neglect of the recoil to the atom during scattering 
is valid p9|. In the first Born approximation we ne- 
glect multiple scattering events and calculate the total 
scattered field as the sum of the fields scattered by in- 
dependent single scatterers out of the probe field. This 
approach is justified as long as the contribution from all 
other scatterers to the local field seen by an individual 
scatterer is small compared to the probe field which de- 
mands e.g. negligible absorption. This condition is also 
helped by destructive interference of scattered waves for 



the case of balanced sub-level populations in the experi- 
mental configuration introduced in the previous section. 
The condition of independent scatterers can only be met 
for not too high atomic density (riat < k^) otherwise 
resonant dipolc-dipole interaction can change the single 
atom scattering properties appreciably. 



A. Scattering integral 

As the first ingredient we need the scattering ampli- 
tude / for a single atom. In our calculations we replace 
the p-wave scattering amplitude by an isotropic ampli- 
tude f equal in magnitude to the true forward scattering 
amplitude. This assumption greatly simplifies analytic 
calculations and will be justified a posteriori by the ob- 
servation that constructive interference occurs only close 
to the forward direction. Integrating |/p over the solid 
angle renders the scattering cross-section at a given de- 
tuning A of the probe laser from the atomic resonance 
|30| . For the experimentally relevant case of alkali atoms 
our choice for the scattering cross section is valid for any 
sub-level of the ground state probed by linearly polarized 
light provided the detuning A is large compared to the 
excited state hyperfine splitting: 
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We remark here that this classical treatment includes the 
full response of the atom to the incident field, i.e. also 
the coupling to all empty modes of the electromagnetic 
field. There is no distinction between spontaneous and 
induced emission and the only assumption made is that 
the response of the atom is linear in the incident field, 
i.e. the atomic transition is not saturated. Naturally, 
since we treat the field as a classical variable, the sta- 
tistical properties of the field, e.g. also the frequency 
spectrum of detected scattered light, are not described 
correctly. However, in the experimentally relevant case 
of large detuning and low intensity scattering is almost 
entirely elastic 31] and in the first Born approximation 
these shortcomings do not enter the problem, since ex- 
change of photons between the atoms, i.e. multiple scat- 
tering, is neglected. 

The total scattered wave is the sum of the waves scat- 
tered by individual scatterers j {j = l..Nat) and can be 
expressed in terms of electromagnetic field vector in com- 
plex notation as [Fig|2]: 

i^..(.o-E4.o.(^-)/ ^^^^:f'^l:^-'^ 
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Figure 2: Huygens-Fresnel propagation used in the model 



The atomic density or density of the point scatterers 
can be written in the form: 

Nat 

n(r) = S{r — -fj) 

with natural normalization condition giving the number 
of atoms in the atomic sample: 



Nat = / dr^n{r) 

In the following we use a continuous density distribu- 
tion according to a smooth probability distribution to 
find a particle in a small volume element, again suitably 
normalized to the total number of particles. This aver- 
aging procedure eliminates large angle scattering off the 
microscopic density fluctuations, which is equivalent to 
the single particle Rayleigh scattering background for a 
discrete distribution of scatterers, i.e. spontaneous emis- 
sion at low saturation parameter Equally Bragg 
scattering of light from a regular distribution of atoms 
on the wavelength scale is lost by the coarse graining. 
The scattered field at some observation point r' outside 
the sample in integral form becomes: 



Esc{r') 



dxdy 



dzEprobe{r)n[r)f p 

L \r' — 7^ 



(6) 

In the paraxial domain (small angles to the optical axis) , 
where we expect constructive interference of scattering 
amplitudes to be concentrated, we can approximate the 
spherical wave propagator in eq.© by using a Fresnel 
expansion formula for the distance |r' — r| 



Since we want to describe free diffraction of probe light 
and scattered light on an equal footing, we choose the in- 
cident probe beam not as a plane wave but rather as gaus- 
sian with parameters w{z), R{z), ^(z) being the beam ra- 
dius, wavefront radius and Guoy phase, respectively: 



Eprobe{x, y, z) = Eq 
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As a realistic model for the density distribution of the 
trapped sample we choose a Gaussian function. This cor- 
responds to the equilibrium shape of a thermal distribu- 
tion of atoms residing inside a harmonic oscillator poten- 
tial. For low enough temperature with respect to the trap 
depth this is a good description for a dipole trapped sam- 
ple. In particular, for the transverse dimensions, where 
we carry out the the integration of the scattering integral 
analytically, this choice simplifies the mathematics. We 
note in passing that a generalization to arbitrary trans- 
verse distributions by expansion into Hermite polynomi- 
als is possible though analytically cumbersome. 

The atomic density distribution in the transverse di- 
rection has a radius Wa, which depends on z due to the 
weaker confinement by the dipole trap laser beam (wave- 
length Xdip) away from its minimal beam waist. In the 
longitudinal direction (along the propagation axis of the 
probe beam) it is described by a I/e-length parameter 
Lo. 
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Finally, the scattered wave field can be evaluated by 
solving the integral 
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in the phase factor, while we use |r' — r | ~ z' — z in the 
less critical denominator. Inserting this we can write the 
propagator in eq.® as: 



irfil -1^ cxp{-ik{z' - z)) ( xx' + yy' 
K[\r' — r\) ~ — ^ exp < IK 
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Here the integration over z is to be taken only over the 
length effectively occupied by the sample, but cannot be 
extended beyond the observation plane. We evaluate the 
scattered field distribution in some distant observation 
plane {M' in Figl^J by carrying out the integration over 
the transverse coordinates of the sample analytically and 
integrating numerically over the length of the sample. 
Using standard software on a desktop PC a scattered 
field profile can be calculated in several seconds allowing 
for fast interactive optimization of parameters. Not sur- 
prisingly for our model assumptions and the choice of the 
density distribution, we find the scattered mode profile 
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to be very close to Gaussian in all of the studied cases 
and we can extract parameters like width and radius of 
curvature by fitting to the corresponding mode profile. 
The scattering efficiency is evaluated by calculating the 
total scattered power in the observation plane. 
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B. Qualitative considerations 



Before presenting results of the numerical calculations 
some qualitative considerations are at hand to train our 
intuition for the results to be expected. First, the total 
scattered power is strictly proportional to the square of 
the number of atoms for a fixed geometry of the sam- 
ple. This is a simple consequence of our continuum ap- 
proximation for the density distribution and may at first 
sight seem disturbing, but is of course natural for co- 
herent scattering, where constructive interference of sin- 
gle scattering amplitudes occurs in phase-matched direc- 
tions. Secondly, far enough away from the sample all 
scattered waves will interfere constructively in the strict 
forward direction, so the on-axis scattered intensity for 
a wide probe beam will be independent of the exact ge- 
ometry of the sample. This implies that scattering effi- 
ciency is determined effectively by the opening angle of 
the scattering cone around the forward direction. Simple 
scaling arguments can be derived for this opening angle 
by looking at Fraunhofer diffraction from transversally 
and longitudinally extended samples. 

Let us first look at the diffraction cone of a short homo- 
geneous sample of width 2wa (FigOJ. To find the angle 
where interference of scattered waves ceases to be con- 
structive we divide the sample in two halves. For a path 
length difference of half a wavelength between the ends of 
a half, destructive interference will occur, giving a limit 
to the opening angle of the scattering cone. In analogy 
to the far field diffraction angle of a Gaussian beam we 
find for a Gaussian source distribution an opening angle 
vtr of: 



Vtr 



tan Vtr 
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(9) 



Narrow samples scatter thus more efficiently than wide 
samples and integrating the angular distribution of scat- 
tered intensity predicts a 1 /w^ dependence on the trans- 
verse width of the sample. 

Next we consider a pencil-shaped atomic sample. For 
this sample Lq ^ Wa- By dividing the atomic sample 
again into two parts (Fig|3J) we can estimate the angle 
at which the longitudinal extent of the cloud causes de- 
structive interference. Introducing the path length dif- 
ference (5, using a small angle approximation and taking 
into account the gaussian apodization we can estimate 



Figure 3: top: Diffraction limited scattering cone for Gaussian 
pancake-shaped sample; bottom: Diffraction limited scatter- 
ing cone for Gaussian pencil-shaped sample 



the opening angle 9l as follows: 



, 1/2 

^ .Trio, 



(10) 



Equating the two expressions for the opening angle 
we can define a characteristic length z^a, the atomic 
Rayleigh range, to compare the influence of the trans- 
verse and the longitudinal extent: 



X 



(11) 



For atomic samples of length Lq comparable or longer 
than Zra, scattered waves from different sections along 
the propagation direction will be mismatched in phase 
and the total scattering cross section will be significantly 
reduced with respect to a short sample with the same 
number of atoms. 

We can construct an approximate expression for the 
total scattered power combining the above arguments in 
order to cast the influence of the experimental param- 
eters sample width, sample length and beam diameter 
into a compact formula. Neglecting for a moment the 
change of transverse spread over the length of the sam- 
ple, the scattered intensity on the optical axis far away 
from the sample where all atoms are phase matched is 
approximately: 



/sc(0,0,2') 



probe 



Sttz'- 



(12) 



as can be verified easily by integration of eq.® in the 
appropriate limit. To find the total scattered power we 
replace the integration over the solid angle by a multipli- 
cation with 7r/2 • d^ff, where the effective opening angle 
Oef f is chosen with the help of the Fraunhofer diffraction 
considerations from above. The extra factor of 1/2 takes 
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a) Wa [|j,m] b) w„ [|im] 



Figure 4: (color online) a) Power of the scattered wave (symbols) vs. the characteristic transverse radius of atomic sample of 
length I/O = 1/im for constant number of atoms and a wide probe beam wo = lOOO/im together with the analytic prediction 
from ea. ll4ll (solid line), b) FWHM of the intensity distribution in the observation plane for the same parameters. 



care of the very close to Gaussian profile of the scattered 
wave. We expect this to be an excellent approximation 
whenever the scattering cone is narrow. In order to model 
the tradeoff between transversal and longitudinal limita- 
tion we have to design a function which takes the value of 
the smaller of the two angles whenever they are grossly 
different. We take: 



3\ 



1/2 



1/2 



with 
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(13) 



Here the transverse limit angle takes into account diffrac- 
tion both due to the sample width as well as due to the 
probe beam width. There is a great deal of freedom in 
the choice of 0e// and different definitions will lead to 
different functional dependencies of the scattering effi- 
ciency on the length of the sample. Our specific choice 
for 0e/ / is motivated by the crossover we observe in our 
numerical calculations for wide probe beams presented 
below. Inserting the above formula we arrive after some 
straightforward algebra at a compact expression for the 
scattered power as: 



- ^ProbeI'at^^3^2. ^ j_ ,„2 /,„2 



1 



m 1 + Wl/W^ ^1 + (Lo/Sra)^ 

(14) 

Here we introduced z^a, the modified atomic Rayleigh 
range by using the definition of 6t from eg 1131 and the 
relation Zra = A/(7r6 
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C. Numerical solutions of the diffraction integral 

Armed with the intuitive arguments and eq. H14|l for the 
scattering efficiency we can now proceed to present some 



of our numerical results. We choose parameters for the 
D2 line of atomic Cs in our calculations (A = 852n'm) and 
keep probe power, detuning and number of atoms fixed 
for all results presented in this section. Fig^ shows the 
total scattered power in the observation plane for short 
samples of varying transverse size. 

The samples are probed by a wide {wq = 1000/im) 
probe beam. The scattering efficiency drops dramatically 
with increasing sample size as expected. A comparison 
with the 1/w^ • l/{w'^ + Wq) dependence from our ana- 
lytical estimate shows perfect agreement. The full width 
at half maximum (FWHM) of the intensity distribution 
in the observation plane (Fig^J)) reflects the interplay of 
source size and diffraction in the propagation of the scat- 
tered wave. In fact, our observation plane is not located 
in the true far field for all source sizes and the observed 
dependence is equivalent to the behavior of the spread 
of varying size Gaussian beams at a fixed finite distance 
from their minimum waist position. 

We investigate now how the scattered power changes 
with the length of the sample. In Fig|SJi we show the re- 
sult for a narrow sample probed by a wide beam together 
with the analytical prediction as a function of sample 
length in units of z^a- The rather good agreement with 
the simple function was our motivation to design the ex- 
pression for 9eff accordingly. In Fig[3l3 we repeat the 
calculation over a larger range of scaled length for various 
transverse sizes of the sample for wide (open squares in 
FigEb) and narrow (open circles in FiglSJi) probe beam. 
The length of atomic sample is scaled here in units of 
Zra and the scattered power is normalized to its value at 
infinitesimally short sample length. 

The estimate with our simple analytical formula is rea- 
sonable also over this larger range showing quantitative 
agreement at the level of 20% for scaled sample length 
up to {Lo/zra) = 8. The fact, that scaled data for a large 
probe beam diameter agree among each other much bet- 
ter than with data from a small beam diameter is under- 
standable from the way ea. (|14|l was derived, i.e. neglect- 
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Figure 5: (color online) a) Scattered power vs. the characteristic length of the atomic sample with atomic waist radius 
Wa = 20^m probed by a beam wq — 1000/xm. Numerical data (symbols) and analytic prediction from eg. 11411 (solid line) are 
shown together, b) Same as in a) for sample width Wa ~ 3, 5, 10, 20^m (squares) and for a narrow probe beam wo = Wa ~ 20^m 
(circles) with the length scaled to Zra- 



ing explicitly the change of the probe beam geometry 
over the length of the sample. To understand the effects 
of changing probe geometry we consider the sample cut 
into thin slices. We can identify each slice as a source 
for a Gaussian beam wavelet, which initially inherits the 
phase profile of the probe beam and develops wavefront 
curvature upon propagation. Scattered wavefronts from 
the back end and the front end of the sample will thus 
have different curvature limiting the overlap to small an- 
gles around the forward direction. Our analytic formula 
works well for a plane probe wavefront, but fails to take 
into account the positive effect of a focused probe beam, 
which imprints wavefront curvature of the right sign to 
enhance wavefront overlap. For this reason the scaled 
scattering efficiency is slightly higher for narrow probe 
beams and intermediate sample length than the model 
predicts. For very elongated samples we see a systematic 
deviation of the analytical model from the numerical re- 
sult and going to extremely elongated geometries we ob- 
serve a different power-law of the decay than suggested by 
the simple analytic model (see Appendix A). We refrain 
here from further tuning of the analytic model, since this 
limiting case is of little interest for coupling all atoms of 
a realistic sample efficiently to the light field and accu- 
rate quantitative data can be obtained easily numerically 
whenever needed, anyway. We define the geometric fac- 
tor as the function describing universally the length 
dependence: 

5L = PsciLo)/PsM ^ . }^ (15) 

The third important parameter which can be varied 
in a real experiment is the probe beam size. A probe 
beam size very much larger than the sample size will 
not be optimum. While the sample is illuminated homo- 
geneously the field strength experienced by the atoms is 



rather low. The dependence with decreasing probe diam- 
eter predicted by eq.l^^ is the result of a subtle interplay 
of increased single atom scattering at higher intensities, 
increased coherent scattering efficiency for the subset of 
atoms which is inside the volume covered by the probe 
beam and decrease of the number of scatterers contribut- 
ing effectively to the scattered field. In FigE^ we show 
the scattered power as a function of probe beam size for 
various sample length. Again we find that the transverse 
probe beam size dependence is described very well by the 
analytic formula for the case of a short sample. For short 
samples it is advantageous to use a probe beam size as 
small as possible. For longer samples both the numerical 
calculation as well as the analytical estimate predict a 
finite probe beam diameter for optimum scattering effi- 
ciency. For reasons already discussed in the context of 
FigEls the analytical estimate starts to deviate from the 
numerical result for probe diameters comparable or less 
than the sample diameter, but works very well already for 
ratios as small as 2. We separate the trivial dependence 
on probe intensity and degree of transverse localization 
(cx Wf^^ • ^) from the observed behavior and define a 
geometric factor qt describing the influence of the ratio 
of beam size to sample size wq / Wa on the sample scat- 
tering efficiency as: 



For the polarization interferometric setup in Fig^ it 
is not possible to adapt the reference wavefront to the 
scattered wavefront separately, so instead of maximizing 
the scattering efficiency only, one needs to choose the 
input beam size, such that the scattered mode has also 
good overlap with the input mode. For short samples 
it is easy to see, that this can be achieved only with a 
probe size much smaller than the sample width, since in 
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Figure 6: (color online) a) Scattered power vs relative size of sample and probe beam for samples {wa = 10/im) of length 
Lo = 1, 400, 738, lOOO^tm (diamonds, squares, stars, triangles) together with the analytic prediction from ea. l|14|l fsolid lines) in 
scaled units, b) Predicted probe beam size to match input and scattered wave as a function of sample length in scaled units. 



this limit the probe traverses a homogeneous region of 
the sample. For longer samples the scattering cone nar- 
rows and one can achieve good overlap also for a probe 
size comparable to the sample width. Equating the far- 
field diffraction angle of the input beam with the effective 
diffraction angle for light scattered off the sample, we can 
derive an expression for the sample length which approx- 
imately matches input and output modes: 



^0 



Wo 



(17) 



The interesting region for the ratio wq /wa is values bigger 
than 1, i.e. probe sizes comparable or bigger than the 
sample size. The predictions of the above equation are 
shown graphically in FiglSja. The values obtained from 
the analytical formula provide good starting values for a 
numerical optimization of this mode-matching problem. 

With our numerical calculations we explored the range 
of validity of a simple analytical estimate for the scatter- 
ing efficiency from samples of different size and found 
quantitative agreement at the 20% level over a large 
range of parameters. 



IV. APPLICATION TO A DIPOLE TRAPPED 
SAMPLE 

The formula and numerical calculation suggest that 
a global optimum for the scattering efficiency exists for 
any sample, which simply consists in placing all scatter- 
ers in one point in space. This optimum is, alas, un- 
physical, because dipole-dipole interaction in that case 
dominates the scattering physics and, more pragmati- 
cally, a trapped sample is subject to density limitations 



because of collision induced heating and losses. In the 
following we study the case of atoms trapped in a sin- 
gle gaussian beam dipole trap. In thermal equilibrium 
the shape of the atomic sample in a single beam dipole 
trap is determined by the focal parameter of the beam 
We take parameters for Cs atoms trapped by 
laser radiation at X^ip = 1030nm at a constant trap 
depth of Uq ~ kB ■ ImK and fixed sample temperature 
of T = IQQ^K. Specifying the dipole trap laser power 
determines then the focal parameter needed to achieve 
the trap depth, and thus also the thermal radius Wa and 
length Lq of the sample. Limiting the peak density to 
npeak = lO^^cm""^ specifies then the number of atoms 
Nat- From Fig|^ we infer that for long samples a probe 
beam size equal to the sample size will be close to opti- 
mum and choose this for the calculation. Restoring all 
prcfactors and choosing a number of incident probe pho- 
tons Uph = 10^ jS^ at a detuning A/T ~ 100 we can 
then numerically determine the achievable SNR accord- 
ing to eq.jSJl at unity quantum efficiency as a function 
of the power of the dipole trap laser. In Fig[7^ we show 
the achievable SNR in this configuration, together with 
the number of trapped atoms. We observe that for big- 
ger samples the SNR approaches a constant value. With 
the constraints we placed on temperature and density, 
the benefit of having more atoms is reduced by the in- 
creasingly unfavorable elongated geometry. It turns out 
that while the aspect ratio of the sample increases with 
increasing size of the dipole trap, the scaled length, i.e. 
the Fresnel number, remains constant at Lg/zra = 14 
js^. With this observation, the higher number of atoms 
is outweighed exactly by the increasing transverse dimen- 
sions of sample and probe beam. 

A realistic model for coherent light-atom coupling effi- 
ciency will have to take into account also the losses due to 
spontaneous emission. In fact, a measurement on atoms 
with a spontaneous emission probability approaching 1 
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Figure 7: (color online) a) achievable coupling strength (filled symbols, left axis) and number of trapped atoms(open symbols, 
right axis) as a function of invested dipole trap power; b) achievable coupling strength (filled symbols, left axis) and probe 
detuning (open symbols, right axis) needed to satisfy 77 = 0.1 (see text for details). 



can hardly be considered non-destructive for the collec- 
tive variable. The result from the 1-D quantum model in 
eq.© predicts actually that the achievable SNR and the 
level of destructiveness are coupled, i.e. the achievable 
SNR is directly proportional to the integrated single atom 
spontaneous emission rate 77 (see also 36, 37J). Since the 
transverse size of the trapped sample changes with the 
invested dipole trap power, rj is not the same for the data 
points in Fig|7|i. In order to make a fair comparison we 
thus calculate for each of the data points the number of 
spontaneously emitted photons per atom by evaluating 
the average intensity experienced by the atoms and inte- 
grating over the pulse time (see Appendix B). In Fig[7|3 
we show the same data rescaled to a probe laser detun- 
ing, such that 77 = 0.1 for every point, together with the 
detuning needed to satisfy the constraint on 77. 

Including the condition of equal level of destructive- 
ness restores the advantage of bigger samples over small 
samples. We note that the absolute numbers for the cou- 
pling strength at fixed 77 cannot be changed by simulta- 
neous variation of the detuning and the number of inci- 
dent photons, since the coupling strength and 77 depend 
in the same way on these two quantities. For the elon- 
gated samples probed by a narrow beam the losses are 
distributed quite unevenly over the sample due to the 
rather inhomogeneous illumination. 

V. RELATION TO EFFECTIVE 1-D MODELS 

At the start of our scattering calculation we expressed 
the coupling strength or achievable SNR in terms of a 
number of coherently scattered photons. This number, 
although convenient to calculate, is not a directly mea- 
surable quantity, since we cannot distinguish coherently 
scattered photons from the incident photons in princi- 
ple and only their interference is observable. The results 
from the previous sections allow us to express this arti- 
ficial scattered power in terms of the incident power and 



the interaction geometry. This also makes a direct com- 
parison to the expression for the coupling strength de- 
rived from the 1-D quantum model possible. Introducing 
the transverse beam area Aph as ttwq and equivalently 
the sample area Aat as ttw? we rewrite ea. (|14() in the 
limit of large detuning as [33 : 

Using eq.(|SJ) we obtain the SNR assuming unit quan- 
tum efficiency detection as: 

Comparing this to the expression obtained from an ef- 
fective 1-D model in section Ull 

^^Natnph(J^^ , (19) 

we see how diffraction effects modify the coupling 
strength with respect to the predictions from the 1-D 
model. The two expressions are similar for the case of a 
sample with effective Fresnel number of 1 and equal probe 
and sample diameter "3^. For our choice of Gaussian 
sample and probe the modification of coupling strength 
can be quantitatively accounted for by the geometrical 
factors. We believe that the asymptotic scaling of the 
geometrical factors will be independent of the specific 
choice of the function describing the shape of the sam- 
ple. We can look at the case of a very elongated sample 
{Lq ^ Zra) probed by a narrow beam of the same size as 
the sample by expanding the geometric factors accord- 
ingly and find in this limit for the coupling strength: 

^3/2 / X \^ / T \ ^ 

^'^ — [2^) "^'-'^"^''(^^j ■ (20) 
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The achievable coupHng strength becomes independent 
of the sample size in this limit and is linear in the peak 
atomic density Upeak, instead of being proportional to 
column density as in the 1-D model . The numerically 
observed scaling with the length of the sample (see the 
appendix) reduces the coupling even more for extremely 
elongated samples. 



VI. CONCLUSION 

In this paper we have outlined an efhcient method to 
include diffraction effects in the coupling of light to col- 
lective variables of atomic samples and applied it to an 
experimentally relevant case of atomic ensembles stored 
in single beam dipole traps. The use of gaussian light 
fields is well adapted to real experimental geometries and 
allows for a largely analytical treatment. Tayloring the 
sample and beam geometry, such that probe mode and 
scattered mode coincide is possible and will be useful for 
polarization interferometry or multi-pass experiments. 

Several approximations have been made, mainly to 
keep the model as transparent as possible, and some of 
them can be lifted in future extensions of the model. The 
leading order effect of multiple scattering and particle 
statistics on the refractive index, which we are effectively 
calculating in a single scattering approximation, can be 
accounted for by a correction term depending on the 
local density of scatterers which will allow to calculate 
the geometry of the scattered mode also for higher 
densities Ull. 



out multiple scattering the change of momentum due to 
scattering introduces correlations between internal and 
external variables ji^l . Also, a focused probe beam with 
inhomogeneous intensity distribution across the sample 
exerts a dipole force on the sample leading to contrac- 
tion or expansion depending on the sub-level populations 
for the setup in Fig^. This leads to an effective de- 
cay mechanism for the macroscopic coherence between 
the sub-levels. Similar effects occur naturally also at the 
level of quantum fluctuations. Ultimately, a proper quan- 
tum model will have to take into account the scatter- 
ing induced dynamics of the density correlation function 
of the sample, which determines the structure factor for 
light scattering |43 . lisj . Prominent examples for the key 
importance of the photon recoil for collective scattering 
are the observation of super-radiant Rayleigh and Raman 
scattering in Bose-Einstein condensates [l^ li^, cavity 
cooli ng [ la. I49I I and collective motion in high-finesse cav- 
ities isol. 
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VII. APPENDIX 



Our model is classical in nature, but the point scat- 
terer model can be used also to analyze quantum noise 
contributions and their dependence on geometry. Giving 
up the continuous density distribution one can determine 
numerically the scattering efficiency from randomly dis- 
tributed samples and this way statistically analyze the 
noise introduced onto the scattered wavefront by density 
fluctuations on different length scales. This models spon- 
taneous emission noise as well as nontrivial effects like the 
inherent mode-matching noise discussed in [l^ . There 
are already studies using a wave function Monte-Carlo 
technique to address the effects of spontaneous emission 
in a quantum description |42j| adapted for our case of 
trapped inhomogeneous samples, but the extremely fast 
increase of the dimensionality of Hilbert space limits the 
treatment to very small numbers of atoms only. A clas- 
sical point scatterer calculation can be used conveniently 
also to model experimental imperfections, e.g. alignment 
errors, where the analytical integration over the trans- 
verse distribution would become much more involved. 

The assumption of infinitely heavy scatterers, i.e. the 
neglect of photon recoil, needs closer attention, when the 
fluctuations of collective variables are studied. In fact, 
when working with collective atomic variables one usu- 
ally assumes that internal and external degrees of free- 
dom of the atomic sample are decoupled. Already with- 



The analytical estimate for the total scattered power 
given in ea. H14|l is seen to fail for very elongated sam- 
ples. This can be traced back to the assumption of ho- 
mogeneous illumination in the calculation of the on-axis 
intensity of the scattered field. A simple way to arrive at 
an improved estimate is to introduce an axial average of 
the incident intensity in order to take into account the 
diffractive spreading of the incident beam over the sample 
length. Since the scattering efficiency depends quadrati- 
cally on atom number, the average is performed over the 
squared density distribution and to simplify the math 
the gaussian atomic density distribution is replaced by a 
rectangular distribution of same peak height and area: 




1 



1 + {Z/Zrf 



-dz 



Zr , ,\/Tr/2Lo 
arctan( — j . 



(21) 



The inhomogeneous axial illumination changes also the 
effective length of the sample entering the estimate for 
the opening angle of the diffraction cone. Incorporating 
this effect we find empirically an improved expression for 
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the longitudinal geometrical factor (see ea. (|15|l ) 



with 



9l = 



■ arctan < 



1/2^ 



(22) 



1 



which fits our numerical data at the level of 20% for Fres- 
nel numbers of the atomic sample up to 80. 



B. 



and 



1 + (^/a)- 



■exp(— = an exp{a^)Erfc{a) 



Within the framework of the point scatterer model the 
distinction between spontaneous and induced emission is 
blurred and with the approximation of a microscopically 
continuous density distribution for the point scatterers 
spontaneous emission is completely lost. Introducing im- 
plicitly microscopic density fluctuations and assuming 
that single atom spontaneous emission happens indepen- 
dently of the presence of the neighboring scatterers the 
spontaneously emitted power can be calculated as: 

{Pspon) = Nata [ I,r,,{r)n{r)dV (23) 

Setting for simplicity the wavelength of the dipole trap 
laser, which determines the change in transverse dimen- 
sions of the atomic sample, equal to the wavelength of the 
incident radiation, the single atom spontaneous emission 
rate can be written in the form: 

_ jPspon) _ ^ 2P be 1 ., ^-1/2 

Y Nat ~~^™0~ i+'^iwa/wa)'^ 

X IZo i+d/aV exp(-z2)dz (24) 



where Erfc{a) denotes the complementary error func- 
tion. The separation of the expression for the coupling 
strength into integrated spontaneous emission i] and 
an effective optical depth a as in the 1-D description can 
be done in principle, but does not lead to simple an- 
alytical expressions. In fact, such a separation is also 
not very meaningful when done globally, since the inte- 
grated spontaneous emission rate can have substantial 
local variations due to the inhomogeneous illumination. 
In addition, the contribution of a central volume element 
to the total scattered field and thus to the scattering ef- 
ficiency is much bigger than for a volume element on the 
rim of the density distribution. This means that a spon- 
taneous emission or optical pumping event in the center 
of the sample leads to a more pronounced change in the 
scattering efficiency. In a quantum description the spatial 
inhomogeneity of both light and atom variables naturally 
suggests importance sampling and leads to a concept of 
collective variables which are no longer fully symmetric 
with respect to exchange of single particle labels [l^ . 
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